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Abstract 

A novel approach to the dynamics of dilute solutions of polymer molecules under 
flow conditions is proposed by applying the rules of mesoscopic nonequilibrium ther- 
modynamics (MNET). The probability density describing the state of the system is 
taken to be a function of the position and velocity of the molecules, and on a local 
vector parameter accounting for its deformation. This function obeys a generalized 
Fokker-Planck equation, obtained by calculating the entropy production of the sys- 
tem, and identifying the corresponding probability currents in terms of generalized 
forces. In simple form, this coarse-grained description allows one to derive hydro- 
dynamic equations where molecular deformation and diffusion effects are coupled. 
A class of non-linear constitutive relations for the pressure tensor are obtained. 
Particular models are considered and compared with experiments. 

Key words: Fokker-Planck equations; FENE models; Polymer solutions; 
Non-Newtonian Phenomena; Constitutive relations. 



1 Introduction 



The statistical study of chain models of polymer molecules [1,2,3] have led to a 
collection of constitutive relations for the description of the elastic properties 
of polymer solutions [4,5], known as FENE dumbbell models, in which the 
restoring force depends on the end-to-end vector of the chain. By avoiding 
phenomena related to the fine structure of the polymer molecules, those models 
depict an intuitive picture of the physics involved in the molecule deformation 
[6], and reproduce non-Newtonian phenomena, such as birefringent pipes [7], 
negative wake [8] and cusp formation [9]. 
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Numerical simulations of molecular chain models have shown that the FENE 
approximation is useful to describe systems under flow conditions, for which 
the largest relaxation time of the polymer dominates the dynamics of the 
molecule, but are inaccurate for molecules under unsteady flow conditions in 
which molecules can explore faster relaxation modes [10]. Although simula- 
tions predict a more realistic restoring force, they are unable to describe its 
evolution in terms of a simple mathematical formulation for arbitrary flows, 
legitimizing the search for a simple description which captures short time scale 
physics. 

In this paper, a theoretical framework originally proposed in Ref. [11,13] has 
been adopted to cope with the main features of dilute polymer suspensions 
under stationary flow. A coarse-grained description of the system dynamics 
is given in terms of a Fokker-Planck equation for the probability density de- 
pending on the parameters accounting for the motion of the polymer center of 
mass, and a single vector parameter (called the distortion vector) accounting 
for its elastic deformation. The Fokker-Planck equation can be constructed 
after invoking probability conservation, calculating the entropy production in 
the phase space of the system and following the scheme of mesoscopic non- 
equilibrium thermodynamics [12]. Moreover, it can be used to formulate a 
hydrodynamic description in terms of the evolution equations for the mo- 
ments of the probability. Application to particular cases can be performed 
by incorporating different models for the Onsager tensor coefficients and the 
finite-size of the molecules [14] . 

The paper is organized as follows. In Sec. 2, the derivation of the generalized 
Fokker-Planck equation is presented. In Sec. 3, the evolution equations for 
the moments of the probability are derived, and in Sec. 4, are applied to 
obtain constitutive relations for the pressure tensor in particular cases. Sec. 5 
is devoted to summarize our main results. 



2 Fokker-Planck equations 

The system under study, a dilute polymer solution, will be regarded as a poly- 
mer 'gas' suspended in a Newtonian solvent under stationary flow conditions 
given by vo(r). The description of this system can be carried out in terms of 
a Fokker-Planck equation. 

As a first approximation, the state of the polymer gas can be characterized 
by the positions r and instantaneous velocities u of the center of mass of the 
polymer molecules [11,13,14,15], and by the vector R characterizing the gross 
distortion and orientation of the molecule [3,6]. Here, it is assumed that the 
molecule segments explore different configurations much faster than the vari- 
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ations in the gross distortion due to the local velocity gradient, i.e. molecules 
are distorted instantaneously with the fluid element in which they are embed- 
ded. Notice that, unlike in the conventional approaches [3], here we introduce 
the velocity u as an independent variable necessary to determine the nonequi- 
librium state of the system at the conditions imposed by the velocity gradient. 
In the dilute regime, the description of the state of the system can be carried 
out in terms of the single molecule probability density /(r, u, R, t), normalized 
to the number of molecules and satisfying the continuity equation 

| + V.(u/) = -A.l/VJ-A.[yv,], (1) 

where and Vr are unknown streaming velocities in u and R subspaces 
[13,14,15]. 

The Fokker-Planck equation describing the dynamics of the system follows 
from the explicit expressions of V„ and \r. These can be obtained using the 
MNET formalism which introduces the Gibbs entropy postulate [16,11,13] 

s{t) ^-ke fc In -^dudR + s^q, (2) 

J Jleq 

where /c^ is the Boltzmann constant, s and Sieq are the entropy per mass unit 
and the entropy at local equilibrium of a given reference state, respectively. 
C = mf / p represents the mass fraction of polymer molecules with a given u 
and R, m the molecule mass and p the mass density field, given by 

p(r, t) = J mfdudR. (3) 



The probability density fieq characterizing the system at the reference state is 
given by 

/.e,(r,u,R) = e^[--(-°)^--^], (4) 



where T is the temperature of the heat bath, pB is the local equilibrium 
chemical potential and |(u — Vq)^ the kinetic energy of diffusion per unit 
mass. The last term in the argument of the exponential function, F • R, is 
the work exerted by the bath on the molecule in order to attain a distortion 
R. Here F is the restitutive force per unit mass of the molecule thought as a 
springy body. 

The Gibbs entropy postulate (2) is consistent with the Gibbs equation [17,13] 
TSs ^Se- pSp~^ - [ pSCdudR, (5) 
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where e is the internal energy per unit mass, p is the pressure and ji the 
nonequihbrium chemical potential. The last term of Eqn. (5) is reminiscent 
of the one corresponding to a mixture in which the components are specified 
by two continuum indexes (u and R). Eqns. (l)-(5) can be used to obtain the 
explicit expressions for the streaming velocities and Vr, by calculating 
the entropy production a of the system in the (r, u, R)-spacc, and taking into 
account the evolution equation for C [15]. The mass continuity equation is 
given by 

l^-pV.v, (6) 



where v(r, t) is the average velocity field of the molecules defined by 

v(r, t) = j uCdndR, (7) 



and the convective derivative given by ^ = ^ + v ■ V. After taking the time 
derivative of the Gibbs equation (5) and integrating by parts, assuming that 
probability fluxes vanish at the boundaries, one obtains the entropy balance 
equation 

ds „ ^ 



In writing this equation, we have identified the entropy fiow 

J, = -kB J{u- v)/(ln/ - l)dudR (9) 

1 f Tfl f 

- - y (u - v)/F • KdudK - — y (u - v)/(u - v)2dudR, 
the entropy production 



771 C 

- ^ y /(u - v) ■ V(u - vo)2rfudR, (10) 

and assumed that the elastic force F is independent of r and u. The incorpo- 
ration of R into the description, introduces elastic contributions to both the 
entropy flow and the entropy production in Eqns. (9) and (10). They repre- 
sent a novel ingredient not considered in previous MNET analysis (see, for 
example, Ref. [13,14,15]). 
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After taking variations of Eqn. (2) and comparing with Eqn. (10), one arrives 
at the exphcit expression of the nonequihbrium chemical potential 

^=M:in/ + l(u-vor + F-R. (11) 
m 2 



Once obtained the entropy production a (Eqn. (10)) and identified the nonequi- 
hbrium chemical potential /x (Eqn. (11)) conjugated to the probability density 
defined over the phase space of the system, one can assume linear relationships 
between fiuxes and forces, of the form [11,13,14,15,16] 



/V„ = -/L„„ • 1^ - fhuR ■ ^ - • (u - vo) • (Vvo)^ , (12) 

fVR = -fhRu ■ ^ - fl^RR ■ ^ - f^Rr ■ (U - Vo) • (Vvq)^ , (13) 



where (Va)^^- = and the Onsager coefficient matrices L^j satisfy the gen- 
erahzed Onsager relations [18]. Notice that unhke the usual framework of 
nonequihbrium thermodynamics given in r-space [17], in the MNET formal- 
ism, the hnear relationships (13) are formulated in the (r, u, R)-space [16,11]. 
This assumption will lead, in general, to nonlinear constitutive relations in the 
physical space r after a contraction of the description [14]. Finally, substituting 
equations (12) and (13) into the continuity equation (1), one obtains 



d_ 
du 

d 



OR 



|.V.,u/,= 

Ci • u - Vo)/ H L„„ • ^ + L„ij • F/ H — 

m on \ m aR^ 

C.-(u-Vo)/+^L^„-|^ + L^«-fF/+^|^'| , (14) 
m an \ m aR; 



which constitutes the Fokker-Planck equation describing the evolution of the 
system at mesoscopic level. We have also defined the tensors 



Ci = Ly„ + • Vvq and C2 = I^ru + l-Rr ■ Vvq. 



(15) 



From these relations it follows that, in general, the transport coefficients 
may depend on the imposed velocity gradient, implying the breaking of the 
fiuctuation-dissipation theorem at mesoscopic level, and incorporating non- 
Newtonian behavior [14,15]. 
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3 Dynamics of the moleculEir distortion under flow conditions. 



The equation describing the dynamics of the distortion of a single molecule 
can be obtained averaging over u the Fokker-Planck equation (14). Assuming 
that the Onsager coefficients do not depend on u, one obtains the equation 



— + V 

dt 



C2 • (vfl - Vo)g + hRR-{Fg + 



ksT dg ' 
m 9R , 



,(16) 



which governs the evolution of g{v, R, t) = / fdu, that represents the proba- 
bility of finding a molecule with a given distortion R at a point r at time t. 
The last term in equation (16) contains the thermal fiuctuations and clastic 
contributions. Here, we have introduced the average velocity = J ufdu. 
The term C2 ■ (vq — vr) represents a drag force on the molecule. We will as- 
sume that it is proportional to R times the characteristic stress given by the 
velocity gradient imposed by the heat bath. Thus, at first order in R, one has 
C2 • {vr — Vo) ~ — C2 ■ R ■ Vvq. Notice that averaging equation (16) over r, 
an equation for the distortion vector distribution resembling known equations 
for dilute suspensions of polymers is obtained [19] . 

Prom equation Eqn. (16), we can compute the evolution equation for the 
average distortion field (R) (r, t) 

{R){r,t) = - [ RgdR, (17) 



and for the second moment field A(r, t) 

A{r,t) = (RR) = - / RRt/ciR, (18) 

Th J 



which can be interpreted as the moment of inertia tensor of the molecule 
[5]. In these equations, n(r, t) = / gdR is the particle density satisfying the 
continuity equation ^ = —V ■ [nv(r,t)], obtained by averaging (16) over R. 
Taking the time derivative of Eqn. (17) and using (16), after integrating by 
parts one obtains 



^ = • [n(A - (R)(R)) • Wo] + C2 • (R) • Vvo - L^.^ • (F), (19) 

where (F) = n^^ ^ YgdR is the average elastic force depending on r and t 
through the average value (R)(r, t). The first term in the right-hand side of 
Eqn. (19) contains the standard deviation (A — (R)(R)) of the R distribution. 
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In the case of a narrow distribution this term can be neglected, reducing our 
equation to the one obtained in Ref. [6]. 

Following a similar procedure, from Eqs. (16) and (18) it is possible to obtain 
the evolution equation for A 



dA 
Itt 



-2[A- Vvo]' = 2^L 



RR 



m 



I - (FR) 



+2 



hRr ■ (Vvo • VvJ 



(20) 



where the super index s stands for the symmetric part of a tensor, I is the 
identity tensor and we have chosen L^^ = I in (15), to obey the material frame 
indifference condition for (20), [19]. Thus, the left-hand side of this equation 
is an Oldroyd time derivative. The first term in the right-hand side represents 
the elastic contribution, comprising the thermal fluctuations of the monomers 



■^n?)' their spatial restrictions due to the fact that they constitute a 
chain ((FR)). The last term in Eqn. (20) is a genuine nonlinear contribution 
in Vvq arising from contracting the description from the mesoscopic (r, u, R)- 
space to the physical r-space. Moreover, in this equation, we have neglected 
contributions arising from higher order moments of the distribution since their 
characteristic relaxation times are at least one order of magnitude smaller than 
the one corresponding to (R) and A (see, for example, [15]). 



3. 1 Evolution equation for the pressure tensor. 



At hydrodynamic level, the presence of polymer molecules modifies the stress 
field of the flow by adding stresses due to elastic and Brownian motion fluc- 
tuations. These contributions can be calculated by averaging our description 
over the instantaneous velocity u and the distortion vector R of the polymer 
molecules and using the definitions 

= m /(u - v) (u - v)fdudK, (21) 



and 

P-^ = y (FR + RF)fdudR, (22) 

where P^ and P^ are the kinetic and elastic contributions of the suspended 
phase to the total pressure tensor, respectively, [20,21]. 
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The equations for the pressure tensor contributions can be obtained from 
the evohition equations for the moments of the probabihty density / and 
the Fokker-Planck equation (14). Contributions from third and higher order 
moments of the hydrodynamic hierarchy have been neglected [15]. For the 
kinetic contribution, using Eqns. (14), (15) and (21), one obtains the evolution 
equation 



K 



dt 



+ 2 • (Ci + Vvo) 



= 2 p(Ci 

m 



Kr ■ VVo)^ 



(23) 



The time derivative in equation (23) accounts for memory effects in a similar 
fashion as Maxwell like models do [19]. For times t ^ [Ci + Vvq]"-^, the 
time derivative in equation (23) can be neglected, giving an algebraic system 
of equations for the components of P"^ that can be solved for any imposed 
flow. Eqn. (23) incorporates contributions to the stresses arising from the flow 
inhomogeneities through Vvq. These contributions have not been considered 
in previous descriptions where Brownian stresses were assumed to be given 
by P'^ ^ -^^pl, [2]. To simplify the computation of the elastic contribution 
to the pressure tensor, in the following we will consider that the elastic force 
F is proportional to R with a proportionality coefficient ^o^, where F is a 
dimensionless function of tr{A) giving the form of the spring law, and .^o is 
the characteristic spring restitution coefficient per unit mass. Therefore, the 
elastic contribution is P^ = —^oFpA and its evolution is related to equation 
(20). 

We will now proceed to express the evolution equations (20) and (23) in di- 
mensionless form. In the absence of an externally imposed flow, equation (20) 
suggests that A scales with kBT/m^Q. Therefore, we will scale length with 



kBT/Am^o, time with the inverse of the shear rate 7 and pressure with 77^7, 
where rjs is the solvent viscosity. Since the Onsager coefficient represents 
the mobility, it will be scaled with the inverse of the characteristic friction 
coefficient per unit mass, say /?, [3]. The coefficient L/j^, will be scaled with 
f3/C,o as it charactizes the coupling between the distortion of the molecule 
and the drag forces, see Eq. (13). Thus, the evolution equation (20) takes the 
dimensionless form 



HA 2 r / ~ ■ 

— -2[A-Vvo]^- = -[l^«-(I-FA 



+2D 



L,i,-(Vvo-Vv^)-A]', (24) 



where D = 7/5/^o is the Deborah number. In similar form, the dimensionless 
evolution equation for P^ is given by 

+ 2 (P^ . Wo)^ = I (^C, - W . Wo)' - 2^ (P^ . C,y , (25) 



dt 



8 



where c = pkBT(3/r]m^o is a measure of the concentration of polymers, D* — 
P'^/^o is a dimensionless parameter that compares viscous dissipation and 
molecule relaxation times. We have scaled the friction tensor Ci with (3. Fi- 
nally, notice that the dimensionless elastic contribution to the pressure tensor 
is 

cF 

¥^ = -—A. (26) 
2D ^ ' 



Eqns. (24), (25) and (26) represent the constitutive relations for a "gas" of de- 
formable particles embedded into a Newtonian heat bath. It is worth noticing 
that in this simple model we have characterized the internal dynamics of the 
particle by using the single vector parameter R, known as the gross distortion 
of the particle. In the following section, we will apply these equations to the 
case of a simple shear flows with different models for the restoring force of the 
particles and for the Onsager coefficients. 



4 Applications 



We will now proceed to consider the constitutive relation expressed in terms of 
the evolution equations (24) and (25) for particular forms of the Onsager coeffi- 
cients, corresponding to simplified models of the suspended polymer molecules 
under shear flow conditions. With these simpliflcations, the constitutive rela- 
tion reduces to models that, at first order in Vvq, resemble the known FENE 
models for polymer solutions. 

The cases studied were a dumbbell, a dumbbell with an isotropic friction coef- 
ficient dependent on its gross distortion, and a dumbbell with a non-isotropic 
friction also dependent on the gross distortion of the molecule. Solutions are 
compared with experiments. 

Dumbbells in a simple shear. In this case, the dimensionless Onsager coef- 
ficients Ci, l^RR and l^ur can be assumed to be isotropic and constant. In 
particular we will assume that Ci = I^rr — I and neglect non-linear terms in 
Vvo by taking L^^ = 0. For spheres, 

W = ^I, (27) 



where a is the radius of the sphere, [14]. Taking the radius a as the char- 
acteristic length scale ^JkBT/Am$,Q of the particle one obtains L^^ — 
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Additionally, we will adopt the FENE spring law [3] 

For times t S> [Ci + Vvojj^-^, the time derivative in equations (24) and (25) can 
be neglected, leading to a set of algebraic equations for P*^ and A. Therefore, 
A is a solution of a non-linear algebraic equation and its components are given 
by: A^^ = L'^{1- Ayy) - 2Ayy, A^y = 2AlyD and A^^ = Ayy. With Ayy the 
real solution of 

SD'^Aly + (3 + L')Ayy - = 0. (29) 



The behavior of A^j. and Ayy as functions of D are common to FENE models 
[19]. For large values of D, A^x and Ayy will tend to and zero respectively. 
The solution to the elastic pressure tensor P-^ depends on A and has no com- 
pact expression, therefore it is not shown here. Similarly, the components of 
the kinetic contribution to the pressure tensor are given by 



and 



P^.-PS, = ^(^ + ^\. (31) 



XX yy 2D* V24 D 



whereas the dimensionless first normal stress difference is defined by 

^1 = -D (P^x - P^y + - P^y) . (32) 



The behavior of A^i as a function of D reproduce the quadratic response of 
polymer solutions at low Deborah numbers observed in FENE models [4,22]. 
We have obtained two contributions to the shear stress: a — —D{P^ + P^). 
The contribution of the elastic part is similar to that of the usual FEN13-P 
models [19] . A new contribution arises from the kinetic part and corresponds 
to that of a 'gas' of Brownian particles [15]. 

Isotropic non-constant friction. In this case, the interaction between the molecule 
and the bath is modeled through a size dependent friction coefficient. We are 
assuming that the average shape of the molecule ranges from a sphere in qui- 
escent liquid to a symmetrical spheroid with a principal axis L in the case of 
maximum distortion. In such case, the friction coefficient Ci will be a tensorial 



10 



quantity. However, the ellipsoid will be considered to be aligned with the flow 
and, under this condition, the friction coeflicient will be taken to be that of a 
spheroid in an axisymmetric unbounded flow given by 

^ (a^ + 1) coth""^ a — a' ^ ^ 



where ISq is now the Stokes friction coefficient, a the principal spheroid axes 
along the flow, b the perpendicular one and a = a/y/a'^ — 6^, [23]. Since in our 
formulation the molecule deformation is solely given by R, then a ~ ||R| — 

l^tr(A). Hence, the maximum value of a is |L and as a consequence, a ranges 
form Ljsjl?- — 4fo2 > 1 in the case of maximum deformation, to infinity in the 
case of a sphere (no deformation). Prom this consideration, it follows that 
coth~^ a is never singular. Then, approximating coth~^ a ~ ^ + equation 

(33) becomes 

/? c ^/^o^MA) 

Notice that, by definition, h must be a decreasing function of tr{A), therefore, 
the second term in the approximation (34) must be negligible small as soon as 
the molecule model begins to experience some distortion. Now, by using Eqn. 

(34) the dimensionless evolution equation for the distortion tensor A takes the 
form 

§-A.Vvo-Vvo^-A=— E=1(I-FA), (35) 
dt 9^tr(A) D ^ ^ 

which constitutes a correction to the friction equivalent to those used in Refs. 
[24,5,6]. Notice that the dependence of the friction coefficient on A introduces 
a non-linear coupling between Eqn. (35) and the corresponding evolution equa- 
tion for P^. Solutions for A and can be obtained by following the procedure 
of the previous subsection. 

Figure (1) shows a fit of the isotropic friction model to simple shear experi- 
ments with dilute micellar solutions in water [25]. Experiments were performed 
increasing the shear rate 7 from zero in steps, the steps in shear rate were taken 
at an acceleration of 0.66s~^. For shear rates over 100s~^ the model overpre- 
dicts stresses. The kinetic contribution is responsible for the agreement for 7 
below 25s~^ whereas the elastic contribution dominates for larger shear rates. 

Non-isotropic friction. Here we will take the friction tensor Ci to be diagonal 
but non-isotropic. This case represents a more realistic model than the pre- 
vious ones because it better accounts for the drag on a body with cylindrical 



1 + 



452 



3tr{A) 



(34) 
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symmetry oriented with the flow, [23] . Thus, in a shear flow on the x direction 
^RR = = -^fCi ■ Ci and — for simplicity Where the friction 

tensor will be modelled by 




1 
1 



(36) 



Figure (2) shows a fit of the non-isotropic friction model to experiments with 
the same system [25] but in this case decreasing the shear rate 7 to zero in 
steps at the same acceleration. For this model elastic stresses dominate at 
all shear rates. The models suggests that the difference between experiments 
arises from the possible configurations that the particles can explore in each 
case. The isotropic friction model best fits the stress behavior in the coil- 
stretch transition that occurs during shear rate growth. The residual stress at 
zero shear rate obtained in experiments (Figure (2)) suggests that particles 
have not reached a coiled configuration. This explains why the non-isotropic 
friction model reproduces the stress behavior in shear rate reduction. 

Notice that, in this description, the dynamics of the molecule segments is not 
included. Therefore, the effect of hydro dynamic interactions between segments 
are not contained explicitly in the set of Onsager coefficients appearing in 
equations (12) and (13). If considered, these interactions will contribute to 
the rhcology of the model by modifying the constitutive relations through 
a different set of Onsager coefficients. The aim of the present model was to 
estimate the collective effect of the molecule segments through a set of effective 
Onsager coefficients [27]. 



5 Conclusions 



A novel approach to the description of the dynamics of dilute polymer so- 
lutions has been explored using the formalism of mesoscopic nonequilibrium 
thermodynamics in the simplified case of a single vector parameter character- 
izing the polymer configuration, common to FENE models. 

The description has been carried out by formulating the laws of nonequilib- 
rium thermodynamics after assuming valid the hypothesis of local equilibrium 

in phase space. By calculating the entropy production of the system, we have 
obtained general constitutive equations at mesoscopic level for the dissipative 
currents in r, R and u spaces. Using then probability conservation, a general- 
ized Fokker- Planck equation describing the dynamics of the system has been 
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obtained. 



A general closure to the governing equations for the hydrodynamics of this 
system model is derived from the Fokker-Planck equation. Those equations 
consist in two evolution equations for elastic and kinetic contributions to 
the pressure tensor. The evolution equation for the elastic stresses resembles 
known FENE models. The kinetic contribution, having no counterpart in the 
FENE models, enters into the hydrodynamic description due to the coupling 
between elastic and Brownian effects in the mesoscopic constitutive relations. 

This method can generate a family of constitutive relations in configuration 
space depending on the modeled Onsager coefficients. Three examples were 
analyzed in the case of simple shear flow. In the flrst one, we have analyzed the 
dumbbell model with constant friction recovering the elastic behavior found 
in the FENE models [4]. In the second case, the friction of the molecule has 
been modeled with that of a spheroid aligned with the flow [24,26]. In this case 
we have obtained a steeper coil-stretch transition and good agreement with 
experiments of shear rate growth. In the last example we have considered the 
tensorial character of the Onsager coefficients by taking the non-isotropic fric- 
tion tensor of a spheroidal particle aligned with the flow. This model is capable 
of describing the shear stress behavior in shear rate reduction experiments. 

Given that the general constitutive equations here obtained can be reduced, 
with simple arguments, to models that predict experimental data, we are con- 
fldent that further analysis can provide simple mathematical expressions for 
constitutive relations better approaching more complex situations. For exam- 
ple, in the case where conflgurational rates of change are relevant, the time 
change in the gross distortion of the molecule ought to be included in the 
description of the system, introducing an additional relaxation time. 
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20 40 60 80 100 

Fig. 1. Fit of experimental data (points) for the non-Newtonian shear stress ver- 
sus the shear rate, using the isotropic friction model (line). The fit is made using 
cr] = OmSSkg/ms, {3 = 3.5s-\ ^ = 11 s''^ and L = 80. 
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Fig. 2. Fit of experimental data (points) for the non-Newtonian shear stress versus 
the shear rate, using the non-isotropic friction model (lines). The fit is made using 
CT] = 0.025kg/ms, (5 = 5s~^, ^ = lls~^ and L = 110, 140, 170 starting from bottom. 
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